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ABSTRACT 

We present results from three Suzaku observations of the z = 3.91 gravitationally lensed broad 
absorption hne quasar APM 08279+5255. We detect strong and broad absorption at rest-frame en- 
ergies of < 2 keV (low-energy) and 7-12 keV (high-energy). The detection of these features confirms 
the results of previous long-exposure (80-90 ks) Chandra and XMM-Newton observations. The low 
and high-energy absorption is detected in both the back-illuminated (BI) and front-illuminated (FI) 
Suzaku XIS spectra (with an F-test significance of >99%). We interpret the low-energy absorption 
as arising from a low-ionization absorber with log A^h^23 and the high-energy absorption as due to 
lines arising from highly ionized (2.75 < log ^ < 4.0; where ^ is the ionization parameter) iron in a 
near-relativistic outflowing wind. Assuming this interpretation we find that the velocities in the out- 
flow range between 0.1c and 0.6c. We constrain the angle between the outflow direction of the X-ray 
absorber and our line of sight to be <36°. We also detect likely variability of the absorption lines (at 
the >99.9% and >98% significance levels in the FI and BI spectra, respectively) with a rest-frame time 
scale of ~1 month. Assuming that the detected high-energy absorption features arise from Fe xxv, 
we estimate that the fraction of the total bolometric energy injected over the quasar's lifetime into 
the intergalactic medium in the form of kinetic energy to be > 10%. 

Subject headings: cosmology: observations — X-rays: galaxies — galaxies: active — quasars: absorp- 
tion lines 



1. INTRODUCTION 

Recent observations of nearby galaxies indicate 
a Mbh-Q- relation (e.g., iFerrarese fc Merritt I l2000t 
iGebhardt et al.ll2000( l. where Mbh is the mass of the 
central black hole and a is the velocity dispersion of the 
stars in the bulge of the host galaxy. The presence of 
a MBH~cr relation suggests that a feedback mechanism 
exists regulating the co-evolution between the massive 
black hole at the center of a galaxy and the forma- 
tion of its bulge. A possible mechanism of feedback is 
quasar outflows. Recent theoretical models demonstrate 
that quasar feedback can serve as a fundamental ingre- 
dient in structure formation and galaxy mergers (e.g., 
Granato et"al] 120041 : iHopkins et alTl2005l : ISpringel et all 



2005() . Quasar outflows could possibly provide an impor- 



tant source of feedback during the growth of the super- 
massive black- holes (SMBHs) in galactic bulges (e.g., 
iFabian I [1999). Another possible mechanism of feed- 
back may be linked to powerful jets apparently driven 
by magnetohydrodynamic forces. As observations in- 
dicate, these powerful jets are predominantly present 
in radio-loud (RL) AGNs,^ which show a tendency to 
be fo und in massive g alaxies and dense environments 
(e.g., iBest et aD l2005| 1i. The importance of jets as a 
feedback mechanism has been demonstrated with recent 
Chandra observations of cavities in clusters of galaxies 
and giant elliptical galaxies (e.g.. iMcNamara fc Nulsenl 
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^ Radio-quiet (RQ) AGN in general do not contain large (i.e. 
kpc) scale coUimated jets, although p c-scale jet s have been found 
in some RQ AGNs (e.g., BlundcU et^aD HQM ). Also a fraction 
(~40%) of radio-quiet AGN could have kpc radi o-structures pos- 
sibly indicating the presence of an "aborted jet" (|Gallimore et afl 
[2009) . 



l2007l and references therein). The injection of power 
into the Intergalactic Medium (IGM) by radio jets is 
a promising feedback mechanism that may explain the 
suppression o f cooling flows in the centers of clusters of 
galaxies (e.g.. IFabian et al.ll2000l: iMcNamara et~aI1 l2000l: 



iSchindler et al.ll2001t iHeinz et al.ll2002D ~ It is not clear, 
however, if radio jets also contribute to the feedback 
process in field galaxies, especially ones in the redshift 
range of z = 1 — 3 where the number density of galaxy 
mergers is thought to peak. Most clusters of galax- 
ies are not for med until z k, 1 as inferred from ob- 
servations (e.g.. iHilton et al.ll2007D and as predicted in 
theories th at consider a low-den s ity (17^ ~ 0-3) Un i- 
verse (e.g., iBahcall fc FmTI [l998l : lYounger et all l2005h . 
In addition, the fraction of radio-lo ud AGNs (RLF) ap- 
pears to evolve with redshift (e.g . . iPeacock et al] 119861 : 
ISchneider et al.lll992t iJiang et al. l2007f) and l uminosity 
(e.g.. iLa Franca et all 11994 iJiang et al.l[2007| ). In par- 
ticular, the RLF tends to increase with luminosity and 
decrease with redshift. For example, for luminous AGNs 
(^'^2500 = —26; where Af2500 is the absolute magnitude 
at rest-frame 2500A) it is expected that the RLF declines 
fr om 24.3% to 4.1% as the redshift increases from 0.5 to 
3 (IJiang et al.ll2007h . ^ 

Quasar outflows present a promising mechanism of 
feedback in high-redshift quasars and possibly in both 
radio-quiet and radio-loud AGNs. Powerful winds are ob- 
served in Broad Absorption Line (BAL) quasars, which 
show deep and broad absorption features from highly 
ionized ultraviolet (UV) transitions. BAL quasars are 
also commonly detected to be X-ray weak as a result of 



3 As in IJiang et all ||2007|) RLF can be written in the form of 
log[RLF/(l - RLF)\ =bQ + 6zlog(l + z) + byi{M25oq + 26), where 
M2500 is the absolute magnitude at rest-frame 2500 A, feo ~ —0.13, 
bz ~ 2.05, and ~ 0.18. 
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high intrinsic absorption column densitie s (-/Vh) typicaUy 
in the range of (1-50) X 10^2 cm-2 fe.g.. lGallagher et all 
I2OO2I 120061 ). We note, however, that a recent survey 
of BAL quasars obtained from t he cross correla t ion o f 
SDSS and 2XMM cathalogs by iGiustini et all ([2008h 
finds no or lower than typical intrinsic X-ray neutral 
absorption from that found in optically selected BAL 
quas ar samples. In the or ientation-based BAL model 
(e.g., iWevmann et al.lll991h quasar winds exist in most 
quasars; however, because of the relatively small open- 
ing angles of these outflows only a fraction of radio-quiet 
quasars have detectable BAL features in their UV and/or 
optical spectra. Models based on numerical simulations 
and observations suggest tha t the winds of BAL q uasars 
are n early equatorial ( e.g., iMurrav et al.l Il995t lElvis I 
[2000t iProga et all l200d ): however, there are a few ob- 
served cases of BAL qua sars with outflowi ng absorbers in 
the polar direction (e.g.. lZhou et al|[2006f) . Recent stud- 
ies indicate that BAL q uasars comprise ^^15-40% of the 
quasar po pulation (e.g., Chartas 2000; Hewett &l Foltz I 
[2003.,; Gibson et al.L 2008: Dai et al. 2008.). 

Our current understanding of AGN physics suggests 
that the most likely mechanisms to explain the origin 
of outflows in AGN are thermal driving, radiation driv- 
ing (line and continuum), and magnetic driving. Ther- 
mal driving will produce slow winds (with speeds simi- 
lar to the sound speed) at large radii (~ 10'*i?s; where 
Rs = 2GM/c^ is the Schwarzchild radius) and with 
a relatively sma ll ma ss-loss rate {/ ^ O.lMoyr"^) (e.g., 
iBegelman etallll983t Ikrolik et al.lll986f ). Therefore, in 
AGNs thermal driving is not expected to produce fast 
and massive winds and consequently it is likely not a 
major contributor to feedback. 

Given the typical low temperatures of AGN accretion 
disks (r ~ 10^ K) and the large gas densities at the base 
of winds we expect that initially the absorbing material 
will have a relatively low ionization parameter. For such 
conditions radiation-driving can lead to significant accel- 
eration of the absorber. Magnetic driving could also be 
present in strong AGN winds, through the action of MHD 
(magnetohydrodynamic) forces (e.g., Everett 2005). In 
general, we expect MHD and radiation-pressure forces 
to act jointly with the contribution of radiation pres- 
sure becoming increas ingly importan t in so urces with 
higher L/LEdd fe.g.. lEverettl l2005l [20071) . Dust in 
the outflow could also boost the radiation pressure de- 
pending on the spectral energy distribution (SED) and 
column density of the material surround ing the AGN 
(|Laor fc Brandt 1 120021 : iFabian et al.ll2008f) . At the mo- 
ment, evidence for the presence of near-relativistic out- 
fiows in AGN is accumulating (e.g., TChartas et al. 2002 
Reeves et a l. 2003; Pounds et al. 2003; Dadina & CappT] 
2004 IChartas et al.. .2007a: .Zheng fc Wang ,, .20081 ^ ; 
however, there is no satisfactory model that can 
produce outflo ws with the near-r e lativistic velocities 
observed (e.g iMurrav et all 119951: iProga et all l2000l : 
Everett I I2005D . We note that recent studies (e.g., 
Chelouche fc Netzer I [20031 lEverett I [20051 ) indicate that 

* A recent paper by Vaughan & Uttley (2008) suggests that 
some of the claimed near-relativistic outflows, especially in cases 
with narrow absorption lines, are detected at moderate significance 
levels and may be spurious. We note, however, that the statisti- 
cal significance of the blushifted broad X-ray absorption features 
detected in APM 08279-1-5255 and PG 1115+080 is not disputed. 



with the appropriate shielding, initial density of the 
wind, AGN SED and L/LEdd? the efficiency of the out- 
flows can be significantly increased and the outflow ve- 
locities may approach near-relativistic values. 

Due to their high intrinsic absorption, many 
BAL quasars appear as faint X-ray sources (e.g., 
iGreen fc Mathur I (l996l: [Gallagher erall 119990 . Partly 
because of this faintness, it is difficult to detect BALs 
in X-ray spectra, and as a consequence, there are only 
a few cases where X-ray BALs have been detected in 
gravitationally lensed BAL quasars where the magni- 
fication effect results in increased signal-to-noise ratio 
spectra. Observations in X-rays of the BAL quasar 
APM 08279-1-5255, the mini-BAL quasar PG 1115+080, 
and perhaps the low-ionization BAL quasar H 1413—117 
have suggested the presence of near-relativistic out- 
flows of X-ray absorbing mate rial in these objects 
([Chartas et al.ll2002l . I2003l l2007al[bl ). The reported vari- 
ability of the high-energy absorption features is over 
rest- frame time-scales of 1.8 weeks in APM 08279+5255 
(significant detection of variability) and 6 days in PG 
1115+080 (marginal detection of variability) . The anal- 
ysis of these high-redshift quasars implied that outflows 
should have a significant impact in shaping the evolution 
of their host galaxies and in regulating the growth of the 
central black hole. These observations are particularly 
important because they allow us to probe quasar winds 
at times close to the peak of the comoving number den- 
sity of luminous quasars. 

In this paper we describe the analysis of three re- 
cent Suzaku observations of the lensed BAL quasar 
APM 08279+5255. A ~100 ks observation of 
APM 08279+5255 was performed starting on 2006 Oc- 
tober 12 (OBSl), a ^100 ks observation was performed 
starting on 2006 November 01 (0BS2), and a -120 ks 
observation was performed starting on 2007 March 24 
(0BS3). 

Unless stated otherwise, throughout this paper we 
use CGS units, the errors listed are at the l-cr level, 
and we adopt a flat A-dominated universe with ~ 
70 km s-i Mpc~\ f^A = 0.7, and VLm = 0.3. 

2. DATA ANALYSIS 

For the reduction and analysis of our observations we 
used the Suzaku software version 7, which is included 
in HEASOFT version 6.4. To analyze data from the 
X-ray Imaging Spectrometer (XIS) and the Hard X-ray 
Detector (HXD) we used calibration files that are part 
of the Suzaku CALDB database released on 2008 April 
01. 5 

2.1. XIS data analysis 

Our data reduction followed the procedures recom- 
mended by the Suzaku team for Spaced-Row Charge In- 
jection (SCI) data. The data reduction was performed on 
the event files of each XIS instrument (XIS 0, 1, 2, and 3), 
and began with recalculating the PI® values of the unfil- 
tered event files using the XISPI routine. Once the event 

^ CALDB version 20080401. 

^ Each event has a measured "Pulse Height Amplitude" (PHA). 
A calculated "PHA Invariant" (PI) value is obtained using the 
PHA in combination with the instrumental calibration and gain 
drift. For the XIS, the PI column name is "PI", which takes values 
from to 4095. The PI vs. energy relationship is the following: 
E[eV] = 3.95 X PI [channel]. 
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TABLE 1 

Log of observations of APM 08279+5255. 



Date 


OBS ID'' 


Telescope 


Instrument 


Exposure 


Net exp 


Net counts'" 


/2-10" 


2002-02-24 


Cha02'* 


Chandra 


ACIS BI 


88.8 ks 




5723±76 


4.3 


2002-04-28 


Has02'' 


XMM-Newton 


EPIC pn 


100.2 ks 




12928±136 


4.0 


2006-10-12 


701057010 


Suzaku 


XIS FI 


102.3 ks 


71.3 ks 


7760±88 


4.2±0.4 


2006-10-12 


701057010 


Suzaku 


XIS BI 


102.3 ks 


71.3 ks 


3046±55 


3.5±0.5 


2006-11-01 


701057020 


Suzaku 


XIS FI 


102.3 ks 


67.9 ks 


7121±84 


3.8±0.3 


2006-11-01 


701057020 


Suzaku 


XIS BI 


102.3 ks 


67.9 ks 


2855±78 


3.5±0.4 


2007-03-24 


701057030 


Suzaku 


XIS FI 


117.1 ks 


86.4 ks 


6059±104 


4.0±0.3 


2007-03-24 


701057030 


Suzaku 


XIS BI 


117.1 ks 


86.4 ks 


3833±88 


3.9±0.3 



^ Throughout this paper we identify the Suzaku observations 701057010 as OBSl, 701057020 as OBS2, and 
701057030 as OBSS.'' These counts are obtained in the 0.6—9 keV observed-frame band and in the 0.4-7 kcV 
observed-frame band for the FI and BI chips, respectively. In each Suzaku observation, J^25% of the FI counts 
and ft:28% of the BI counts are background.*^ The fluxes (in units of 10~^'^ergs cm~^ ^~^) iu the 2—10 keV 
observed-frame band are obtained using the best-fit absorbed power-law model (model 2; §3) in our Suzaku 
observations. The fluxes measured in the BI chips arc on average less than those in the FI chips. This is 
due to a higher half-power-diameter (HPD) of the XISl instrument, compared to the HPDs of the XISO, 
XIS2, and XIS3 instruments.'^ In this Table we identify as Cha02 the 88.8 ks obse rvation of APM 08279-1-5255 
performed with Chandra in 2002 and analyzed in detail in IChartas et all (120021) . We also identify as Has02 
the 100.2 ks observat ion of APM 08279-1-5255 performed with XMM-Newton in 2002 and analyzed in detail in 
IHasinger et al.l I I2002I) . 



files were reprocessed, we used the XSELECT software 
to apply the standard screening criteria (see the Suzaku 
ABC guide^) and obtain "cleaned" event files. The data- 
screening criteria include selecting events corresponding 
to ASCA grades 0, 2, 3, 4, and 6, Earth elevation angles 
greater than 5° (ELV>5), Earth day-time elevation an- 
gles greater than 20° (DYE_ELV>20), exclusion of pas- 
sages through or close to the South Atlantic Anomaly 
(SAA), and cut-off rigidity criteria of >6 ^ (C0R>6). 
As a final step in screening the data we removed hot- 
flickering pixels through the use of the SISCLEAN rou- 
tine in XSELECT. The total exposure time of each XIS 
chip decreased by w20% after the above screening crite- 
ria were applied. Using the clean event files we extracted 
events in a circular region centered on the source with a 
radius of 150 pixels (2.5'). Background events were ex- 
tracted in an annulus centered on the source with an 
inner radius of ^^230 pixels (3.8') and an outer radius 
of ^^430 pixels (7.1'). Our selected background region 
excludes APM 082794-5255 and the calibration sources 
located near the corners of the CCDs. The response 
matrix files (RMFs) and ancillary response files (ARFs) 
were generated using the XISRMFGEN and XISSIMAR- 
FGEN routines which include the correction for the hy- 
drocarbon contamination^ on the optical blocking filter. 

For the front-illuminated (FI) XIS chips (XIS 0, 2, 3) 
we considered events with energies lying in the range 0.6- 
10 keV, while for the back-illuminated (BI) XIS 1 chip we 
considered events with energies lying in the range 0.4-8 
keV. Due to calibration uncertainties near the CCD Si K 
absorption edge at 1.84 keV, events with energies lying in 
the range 1.7-1.95 keV were ignored in the analysis of all 
four XIS chips. In order to assess systematic uncertain- 
ties in the response files, we fitted the Ni Ka (7.470 kcV) 
calibration line of each instrument. We found similar 
positive shifts in the inferred energies of the calibration 

http:/ /heasarc. gsfc.nasa.gov/docs/suzaku/analysis/abc/ 
8 The XISSIMARFGEN routine incorporates the XISCONTAM- 
ICALC routine which is used to correct the observation for the XIS 
optical blocking filter (OBF) contamination. The absorption due 
to these contaminants depends on the X-ray energy, time, detector 
ID and location on the detector. 



lines of each XIS chip ranging from 10 to 20 eV. These 
shifts in energy were not large enough to cause any sig- 
nificant impact on our analysis, and therefore we did not 
attempt to correct them. The net source count rate for 
each XIS chip and each observation was «0.04 counts 
s~^, with a background of !=a30 % of the source rate. 
The spectra obtained on the FI chips were combined 
using the routine ADDSPEC (in HEASOFT FTOOLS) 
in order to increase their signal-to-noise ratio. In Ta- 
ble [T] we have included information relevant to the XIS 
data analysis. Specifically, this table lists the observation 
ID, exposure time, net exposure time (after the screen- 
ing process), net counts (for the FI and BI chips) and 
flux in the 2-10 keV observed-frame (for the FI and BI 
chips) using the best-fitted absorbed power-law model 
(model 2; §3). We also have included in Table [1] in- 
formation from two previous deep X-ray observations of 
APM 082794-5255. These observations correspond to an 
88.8 ks Chandra exposure (seelChartas et al. 2 002) and 
to a 1 00.2 ks XMM-Newton exposure (see lHasinger et alj 
I2002t l. The counts collected by the XIS FI chips for each 
of our observations are comparable to those obtained in 
the Chandra observation. 

2.2. HXD data Analysis 

Similarly to the XIS case the clean event files were 
obtained from the unfiltered event files following the in- 
structions in the Suzaku ABC guide. The screening cri- 
teria are similar to those applied to the XIS instruments, 
specifically, we used ELV>5°, DYE_ELV>20°, exclusion 
of passages close to the SAA, and C0R>6 (units of 
[^]). The HXD-PIN spectrum was extracted from 
the cleaned events file described above. We extracted 
the source spectra from the clean files XSELECT. In or- 
der to estimate non X-ray background (NXB) events, we 
used version 2 of a time-dependent instrumental back- 
ground event file (referred to as the PIN background 
event file) which was provided by the Suzaku team. The 
PIN background event file was generated with a count 
rate that is ten times larger than the real instrumental 
PIN background. Therefore, we increased the effective 
exposure time of our observed PIN background spectra 
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TABLE 2 

Log of PIN HXD Suzaku observations of APM 08279+5255. 



Epoch 


Net exposure 


10-40 keV count rate (10" 


"^cts s-i) 






source NXB 


CXB" 


OBSl 


88 ks 


49.75±0.24 48.17±0.07 


2.30±0.05 


OBS2 


89 ks 


50.37±0.24 46.05±0.07 


2.39±0.05 


OBS3 


103 ks 


47.07±0.21 44.12±0.06 


2.29±0.05 



^ The CXB counts have been estimated from a fake spectrum gener- 
ated using the FAKEIT command of XSPEC with the model given in 
equation (fTl . 

by a factor of ten. The exposure time was corrected for 
dead time using the HXDDTCOR task, leaving an ef- 
fective exposure time of ^90% of the original exposure 
time. The effective exposure time of each observation, 
together with the count rates (10-40 keV) of the source 
and NXB are presented in Table H The NXB does not 
include the contribution from cosmic X-ray background 
(CXB). Therefore the CXB counts (see Table have 
been estimated from a fake spectrum generated using the 
FAK EIT comman d of XSPEC with the following model 
('e.g.. lBoldtl[l987l ): 



CXB{E) 
9.412 X 10- 



= p 40kcV 



/ E 



V IkeV 



-1.29 



cm s sr keV 



(1) 



The response file used to fit the PIN spectra was ob- 
tained from the Suzaku CALDB calibration files. The 
HXD spectral analysis was performed in the 10-40 keV 
energy range. 

3. SPECTRAL ANALYSIS 

In this section we fit the Suzaku spectra of 
APM 08279+5255 with a variety of models using the 
software tool XSPEC version 12. We also fit the spectra 
with more realistic models based on the photoionization 
code XSTAR. In all spectral models we assume a Galactic 
column density of 4.1 x lO-^^cm^^ (Kalberla et al. 2005). 
Most of this section concentrates on the analysis of the 
XIS spectra; however, in the last paragraph we present 
results from the spectral analysis of the PIN spectra of 
APM 08279+5255. 

3.1. XIS spectral fits. 
3.1.1. XSPEC spectral fits. 

Each observation of APM 08279+5255 provides spec- 
tra obtained with the single BI chip (XISl) and the FI 
chips (XIS 0, 2, and 3). Since the responses of the FI 
chips are similar we co-added the FI spectra from each 
observation. We note that there is no XIS 2 spectrum 
of APM 08279+5255 for our third epoch (0BS3) due to 
the failure of the XIS2 chip.^ To fit the spectra using 
statistics we grouped each XIS spectrum with a suffi- 
cient number of counts. The minimum number of counts 
per bin was also chosen to maximize the signal-to-noise 
ratio in each bin without losing the features in the spec- 
tra and to keep a similar number of spectral bins in each 

8 On 2006 November 9, about 2/3 of the 
imaging area of XIS2 became suddenly un usable 
| |http://heas arc. gsfc.nasa.gov/docs/suzaku/news/xis2. html [) . 



spectrum (Ri70). The minimum number of counts per bin 
chosen for the BI chip was 40 for epochs OBSl and 0BS2 
and 50 for epoch 0BS3. The grouping for the FI chips 
was 100 counts per bin for epochs OBSl and 0BS2 and 
80 counts per bin for epoch 0BS3. Note that for epoch 
0BS3 we have increased the binning of the BI spectra 
due to the longer exposure and decreased the binning of 
the FI spectra to compensate for the loss of XIS2. 

We fit the spectra of APM 08279+5255 with the 
following models: 1) power- law (PL; XSPEC model 
wabs*pow) , 2) absorbed power-law (APL; XSPEC 
model wabs*zwabs*pow), 3) ionized-absorbed power- 
law (lAPL; XSPEC model wabs*absori*pow), 4) par- 
tially covered absorbed power-law (PAPL; XSPEC model 
wabs*zpcfabs*pow), 5) absorbed power-law with a notch 
(APL+No; XSPEC model wabs*zwabs*notch*pow), 
6) absorbed power-law with an absorption edge 
(APL+Ed; XSPEC model wabs*zwabs*zedge*pow), 
and 7) absorbed power-law with two absorption lines 
(APL+2ALi°; XSPEC model wabs * zwabs * [pow + 
zgauss + zgauss]). 

The results of the FI and BI fits with the models de- 
scribed are listed in Table [H The error bars of the fitted 
parameters are given at the 68% level (Ax^ — 1). For 
models 2 to 7 w e assume an intrins ic absorber with a 
redshift of 3.91 pownes et all 119990 . The fits using a 
power-law model (model 1) are not acceptable in a statis- 
tical sense. We next fit the spectra of APM 08279+5255 
with the absorbed power-law model (model 2; Table [3]) 
assuming an intrinsic absorber. The _F-test indicates 
that fits with model 2 result in a significant improve- 
ment at the >99% and >99.9% confidence levels in the 
FI and BI spectra, respectively, compared to fits us- 
ing model 1. Fits with model 2 indicate significant 
intrinsic absorption in APM 08279+5255 with a col- 
umn density of log A'h~23. We also fit the spectra of 
APM 08279+5255 with more complex models that in- 
cluded an ionized and partially covered absorber (models 
3 and 4), however, these fits did not result in a signifi- 
cant improvement (i^-test significance <95%) compared 
to the simpler model 2. 

Fits to the spectra of APM 08279+5255 with mod- 
els (models 5-7) that account for the absorption found 
between 7-12 keV in the rest-frame result in signifi- 
cant improvements (the i^-test indicates improvements 
at >99.9% and >99% confidence in the FI and BI spec- 
tra, respectively) compared to fits with models that do 
not include this high-energy absorption. We note that 
the absorption feature at 7-12 keV in the rest-frame cor- 
responds to a signi ficant detection following t he criteria 
described in §3 of IVaughan fc Uttlev I (|2008f l. Specifi- 
cally, we find the ratio of the total equivalent width^-'^ 
of the absorption features to their uncertainty to be 
EW/aEw ^ 3 in every observation (see models 5-7 in 
Table [31). 

We note that if we replace the APL+2AL model by the 
XSPEC absorption-line multiplicative model wabs * zwabs * gabs * 
gabs (pow), we obtain similar results for the fitted energies and 
equivalent widths of the absorption features found at energies 
7-12 keV in the rest-frame. All the results described in this paper 
using the APL+2AL can be reproduced using this multiplicative 
model. 

" The equivalent width (EW) is defined as EW = J ^""^^^ dE, 
where Ec is the continuum flux and Eg) is the flux in the absorber. 



TABLE 3 

Results from spectral fits to the three Suzaku observations of APM 08279+5255. 



Model" 


Parameter 


Values OBS 1 


FI spectrum'' 
Values OBS 2 


Values OBS 3 


Values OBS 1 


BI spectrum'' 
Values OBS 2 


Values OBS 3 


1 


r 

v2 /i> 


, 70+0.02 
^- ''-'_0.02 
118 3/73 

6 X 10-* 


1 7C+0.03 
^- '^-0.03 
1 2"! 8/66 

-LtrtJ.O / \J\J 

1 X 10-5 


1 77+0.03 
'-0.03 

99.1/71 
0.02 


1 -57+0.04 
114.0/71 

9 X 10-4 


J 50+0.05 
140 S/6fi 

±.'-t\J.Kj/ \J\J 

2 X 10-''' 


1 (51+0.04 

mo 6/71 
0.01 


2 


r 

log A^H 


09 OO+0.09 

90.5/72 
0.07 


1 Qo+0.04 
-'-■y°-0.04 
22 92+"-"-' 

ZZ.yZ_g j^g 

87.9/65 
0.03 


1 SS+O-05 
-'-■S°-0.05 

22.66l°:i5 
91.0/70 
0.05 


1 q6+°°'^ 
^^•'^0-0.07 

92 74+0.09 

'1-0.09 

86.7/70 
0.59 


-1 93+0.07 
l-'''-'-0.07 
99 vo+0.09 

101.3/65 
3 X 10-3 


1 Q9+0.06 

i-'^^-o.oe 

99 7r-+0.09 

ZZ. /0_g gg 
59.0/70 

0.82 


3 


r 

log A^H 
log I 


1 go+^-o" 

99 Q9't"t).08 
^^■"^-0.07 

<0.18 
86.0/71 
0.11 


1 99+0-05 
-^■="^-0.05 
99 99+0-07 

ZZ.3»_g gg 

<0.11 
83.2/64 
0.05 


1 90+0-05 
22.80l°:13 
<0.75 
87.6/69 
0.06 


1 97+0.07 
-0.07 

22.79t°-}i 
<0.22 
84.4/70 
0.12 


1 94+0-08 

-^-■^^-0.08 

22.87l°:i° 
<0.68 
99.7/64 
3 X 10-3 


1 94+°-°^ 
22.93lj',-"!? 
<0.98 
56.6/69 
0.86 


4 


r 

log iVn 
CF 


1 90+°"'^ 
-^■^^-0.07 
99 ocr+0.18 
ZZ.OO_Q j^O 

86.9/71 
0.10 


2 00+° °'^ 
^■'^'-'-0.08 

23.0lt°:f, 
0.89«;^^ 
85.6/64 
0.04 


1 91+0.09 
-^■''-^-0.09 
22.7ll°:22 

92+0-08 
"^■"^-0.26 

88.9/69 
0.05 


1 96+°°* 
-^•^"-0.09 

22.75t°:l^, 
0.951°:°? 
85.4/70 
0.10 


1 02+0.09 
l-='^-0.08 

99 79+0.12 
'=^-0.10 
94+0-06 
'^■='^-0.08 

100.9/64 

2 X 10-3 


1 91+0.06 
l-''l-0.06 
99 7c; + 0.12 
^^•'^-0.07 
93+0-07 
^■="•^-0.08 

57.9/69 
0.82 


5 


V 

log A^H 
-Bnotch[keV] 
Wnotch[keV] 
fnotch[keV] 
EW„oteh[keV]'= 

P(xVu) 


99 nK+a.oa 

9 sg+'^-i'^ 
^•■J^-o.is 

3.98+°-i 

0.80 ±0.21 
62.3/69 
0.70 


1 QS + O-04 
-'-■y°-0.04 

22.98l°;°| 

9.75i;;-i^ 

3-86±'^15 
n 10+0.05 

^■J^°-0.05 

0.70 ± 0.22 
70.6/62 
0.21 


1 qi+0.05 
l-yi-0.03 
99 09+0.10 

9.4U_o 2g 

4-03+r39 

Q 17+0.05 

0.69 ±0.21 

74.2/67 
0.26 


1 95+0. 06 

22.80+|j-E!? 
9 SI +0-26 
y-°l-0.25 
4 7-1+0.88 

26+°'°'5 

'J.^D_g gg 
1.22 ± 0.33 

58.1/67 
0.77 


1 91+0.07 
l--'l-0.07 
99 OQ+0.08 

ZZ.OO_0 gg 

10.411":;;° 

91+0 05 

1.05 ±0.36 
90.8/62 
0.01 


i-92lS;;j^ 

22.801°:°^ 
9.96l°g-3^ 
4.52l°ii 

r, 10+0. 05 

^•J^°-0.05 

0.81 ±0.28 

48.1/67 
0.96 


6 


r 

log ATh 

-Ezedge[keVl 

"^zedge 
EW.edge'^ 

Pix'/u) 


-^•^•^-0.04 

99 94+0.07 

ZZ.»'±_g gg 

7 71+0.15 
'■'-^-0.14 
f, 09+0.07 

1.19 ±0.28 
67.1/70 

0.58 


2 00+° °'^ 
^■'^'^-0.04 

23 nn+o-oo 

Za.UU_g gg 

8 06+° " 

i3.UD_Q 29 
,(,+0.08 

1.12 ±0.28 
71.1/63 

0.23 


1 93+0.05 
^■^•^-0.05 
22.82t°:ll 
740+0.12 
'•^°-0.13 
01 +0.08 

1.09 ±0.26 
74.3/68 

0.28 


1 94+0.07 

-l-^^-0.07 
22.80l°:°« 
7 47+0.13 

'■^'-0.13 
56+0-18 

1.86 ± 0.59 
64.3/68 
0.60 


1 R«+0-07 

l-*8_o gy 

99 SI +0-08 

ZZ.01_g gg 
Q 09 + 0.52 

°"^^-0.33 

46+°-" 
1.54 ± 0.64 
91.0/63 
0.01 


1 90+°-"° 

l--'"-0.06 
22 79+0.07 

ZZ. I»_g gg 

7 74+0.20 

' • '^-0.23 
41+° " 

1.42 ±0.42 
46.1/68 
0.98 


7 


r 

log AfH 
•Babsi[keV] 
f^absi [keV] 
EWabsi[keV]'= 

Sabs2[kcV] 
0'abs2 [kcV] 

EWab.2[keV] = 

x^/y 

P{x^/u) 


1 94+0.04 
-^•^^-0.04 
99 94+0-08 

ZZ. J1_g gj^ 

8.ioi;;:}j, 

0.291°-? 

r, 99+0.08 
"•^^-0.06 
1 60+°'^* 

o.78i;-^ 

49+0-2.5 
^•^="-0.24 

58.0/66 
0.75 


2 02+"-°'^ 

^■^^-0.03 

22 99+»-°^ 

ZZ...IJ_g g4 

8.62l»;?° 
07+0.20 

^■■"-0.18 

0.25t»:i° 
10.851°:- 

n oc+0.21 

97+0-16 
'-'■^'-0.18 

66.4/59 
0.24 


i-94i;;:°^ 

22.8ll°:i° 

7 o-i+O.lO 
' ■"'^-0.09 
<0.32 

21+°-°-' 
'^■^-1-0. 11 

io.26l°:;{^ 

n 72+0.52 
'^-0.33 

46+°-" 
68.8/64 
0.32 


1 95+0.07 
1-^^-0.06 
99 0-1 +0.08 

ZZ.O±_Q gg 

7.861°:?° 

<0.48 

23+°-°'^ 
"•^■'-o.io 

10 53+0-29 

iU.OO_g gg 

97+0-54 
^^•^'-0.52 
-1 19+0.36 
1-1^-0.38 

54.0/64 
0.81 


I 91+o.Ofi 

22 79+0-10 
'-^-0.10 

9 26+°-^'* 
''•^"-0.26 

<0.52 

34+0-22 

'^■■'*-0.20 

II 32+°°-^ 

ll.,DZ_g gg 

<0.92 
<0.98 
79.8/59 
0.04 


l-92l°:°^e 
22.801°:'- 
Q 97+0.32 
"■^ -0.36 

<0.74 
0.39l°:i? 
10.98l°:« 

o.7ii;;:f6 

0.58l°:f, 
43.8/64 
0.97 



^ Model 1 is a power-law with Galactic absorption (PL; XSPEC model wabs*pow); Model 2 is a power-law with Galactic absorption and intrinsic absorption 
(APL; XSPEC model wabs*zwabs*pow) : Model 3 is a power-law with Galactic absorption and ionized- absorption (lAPL; XSPEC model wabs*absori*pow} ; Model 
4 is power-law with Galactic absorption and partially covered absorption (PAPL; XSPEC model wabs*zpcfabs*pow) ; Model 5 is a power-law with Galactic 
absorption, intrinsic absorption, and a notch absorber (APL + No; XSPEC:* model wabs *zwabs *notch*pow) ; Model 6 is a power-law with Galactic absorption, 
intrinsic absorption, and an absorption edge (APL + Ed; XSPEC model wabK*zwal)s*zedge*pow) ; Model 7 is a power-law with Galactic absorption, intrinsic 
absorption, and two absorption lines (APL-1-2AL; XSPEC model wabs * zwabs * [pow + zgauss + zgauss]).^ The spectra fitted are the added (ftools ADDSPEC) 
spectra of the FI chips (XISO, XIS2 and XIS3). For OBS 2 the XIS2 CCD was not operational, and only the XISO and XIS3 spectra were added for this observation. 
The BI spectra are taken with the XISl chip.'' EW stands for equivalent width, which is defined as EW = / — ^-p — —dE, where is the continuum flux and 
is the flux in the absorber. 
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Fig. 1. — Suzaku FI (left panel) and BI (right panel) spectra of the combined images of APM 08279+5255 for the three observations 
(OBS 1, 2, and 3), fit with Galactic absorption and a power-law model to events with observed-frame energies above 3.6 keV and then 
extrapolated to lower energies. In the lower panel of each observation, we show the residuals of the fit with l-cr error bars. High-energy 
absorption features are detected within the 7-12 keV range (dotted lines). We have marked with an arrow the best-fitted energies of the 
first absorption feature of model 7 for epochs OBS2 and OBS3. For the FI spectra the grouping in OBSl and OBS2 is 100 counts per bin, 
and in OBS 3 it is 80 counts per bin. For the BI spectra the grouping in OBSl and OBS2 is 40 counts per bin, and in OBS3 it is 50 counts 
per bin. 



To illustrate the presence of the high-energy absorption 
features, we fit the spectra from observed-frame energies 
of 3.6-10 keV with a power-law model and extrapolated 
this model to the energy ranges not fit (see Figure [1]) . 
The lower panels in Figure [T] show the residuals (differ- 
ence between the measured counts and model) between 
the best-fit power-law model and the FI and BI data, 
respectively. The best-fit values of the photon indices 
in all observations with this model were consistent with 
r = 2 at the l-a level. For the purpose of comparing 
the absorption residuals between epochs the photon in- 
dices for all observations were set to F = 2.0. From these 
fits we notice that the residuals show an absorption fea- 
ture centered near a rest-frame energy of ~8 keV and a 
possible second absorption feature near a rest-frame en- 
ergy of ~10 keV. We fit the high-energy absorption fea- 
tures with the models listed in Table [Sj From these fits 
we found that adding to the APL model an absorption 
edge (APL-I-Ed) or two absorption hues (APL-I-2AL) im- 
proves the fits at the >99% confidence level in the two 
sets of spectra (FI and BI) and in each observation. The 
i^-test indicates that we cannot distinguish between the 
(APL-I-Ed) and (APL-I-2AL) models for fits performed 
to the spectra of APM 08279-^5255 in epochs 0BS2 and 
0BS3, since both models fit equally well the 7-12 keV 
rest-frame absorption during these epochs. However, 
fits to the FI and BI spectra of epoch OBSl using the 
APL-I-2AL model provide a significant improvement at 
the 95% and 98% confidence levels, respectively, com- 
pared to fits that use the APL-t-Ed model. It is im- 
portant to note that the APL+2AL model was clearly 
favored over the APL +Ed model in a pre vious 88.8 ks 
Chandra observation (jChartas et al.l 120021 ). Fits to the 
spectra in epoch OBSl with a model that includes an ab- 
sorption notch (see Table [3]) also provide a significant im- 
provement compared to ones using the APL+Ed model. 
These F-test improvements are at the ~98% and ^^99% 
levels of significance in the FI and BI spectra, respec- 
tively. We note, however, that when we compare the 
quality of the spectral fits that use the APL-I-Ed model 
with fits that use either the APL-f 2AL or APL+notch 
models, the F-test may not be a reliable tool. The rea- 
son for the non-reliability of the F-test is that we are not 
comparing nested models (see Protassov et al. 2002 for 
details). In order to check the reliability of the F-test 
for these cases, we performed Monte Carlo simulations 
of 10,000 fake spectra (using the FAKEIT command of 
XSPEC) assuming an APL-j-Ed model. In these simu- 
lations, the energy and optical depth of the absorption 
edge are assumed to be normally distributed around their 
fitted values for epoch OBSl (see model 6 of Table [3]), 
with a standard deviation given by the error bars of the 
fits. All other parameters of the APL-f-Ed model were 
set to their best-fitted values (epoch OBSl and model 6 
of Table [3]). The results of our Monte Carlo analysis are 
presented in Table ID In each simulation we have fitted 
the data with the null model (APL-|-Ed) and the alterna- 
tive model (cither APL-I-2AL or APL+notch). We then 

We only find marginal improvements in fits to the spectra of 
APM 08279+5255 taken in epochs OBS2 and OBS3 using model 7 
(APL-I-2AL) compared to fits using model 6 (APL+Ed). Specifi- 
cally, in epoch OBS2 these improvements are at the 61% and 91% 
significance levels in the FI and BI, and in epoch OBS3 they are 
at the 68% and 50% significance levels. 
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7.5 8.0 8.5 9.0 9.5 

Rest-Frame Energy (keV) 

Fig. 2. — 68% and 95% confidence contours of absorption-line 
fluxes versus absorption-line energies of the first modeled absorp- 
tion line (model 7, Table [Sjl. The upper and lower panel contours 
correspond to fits performed to the FI and BI spectra, respectively. 
Solid and dotted lines are contours for epochs OBS2 and OBS3, 
respectively. 

calculated the value of the F-statistic between the null 
model and the alternative model. In Table [3] we show 
that the p-value, which represents the fraction of simu- 
lated cases with values of the F-statistic higher than the 
actual value obtained from our real data, is similar to the 
null probability of the F-test. We therefore conclude that 
our F-test values are reliable and are approximately rep- 
resentative of the improvement of the alternative model 
(either APL-f 2AL or APL+notch) with respect to the 
nuU model (APL+Ed). 

The results of the spectral fits shown in Table [3] indi- 
cate a change (greater than 1-a ) of the energies of the 
best-fit values of the first absorption line (absl; model 7), 
and in the absorption-edge energy (FEdgo! model 6) be- 
tween epochs 0BS2 and 0BS3 in both the FI and BI 
spectra. This change is also suggested by the residuals in 
Figure [TJ where we have marked with an arrow the best- 
fitted energies of the first absorption feature of model 7 
for epochs 0BS2 and 0BS3. The shift in the energy of 
the first absorption line is indicative of possible variabil- 
ity of the outflow. This change can be seen more clearly 
in Figure [2] where we show the confidence contours of 
the best-fit energies of the first absorption line (APL + 
2AL model) versus its line-flux normalization in epoch 
0BS2 (solid line) and in epoch 0BS3 (dotted line).^^ 

" We note that for the FI spectra of APM 08279+5255 during 
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TABLE 4 

Estimates of the improvement of fits to the spectra of 
APM 08279+5255 using alternative models to the APL+Edge model. 



Alternative model 


Spectrum" 


_F-statistic/null probability'' 








APL+2AL 


FI 


2.59 / 4.4 X IQ-^ 


5.2 


X 10^ 


-2 


APL+2AL 


BI 


3.05 / 2.3 X 10-2 


3.1 


X 10" 


-2 


APL+notch 


FI 


5.32 / 2.4 X 10-2 


6.8 


X 10" 


-2 


APL+notch 


BI 


7.15 / 0.9 X 10-2 


3.9 


X 10" 


-2 



^ The XIS FI and BI spectra of APM 08279+5255 considered in this comparison arc 
taken from epoch OBSl.^ The value on the left of the slash is the F'-statistic and 

IS given by F = "1^,""! The value on the right of the slash represents the 

probability of exceeding the F-statistic based on the F-test.^ The p-value represents 
the probability of exceeding the F-statistic based on our Monte Carlo simulations. 
The value of the F-statistic is obtained by comparing the null model (APL+Edge) 
with the alternative model listed in the first column for fits performed to the spectra 
of APM 08279+5255 obtained in epoch OBSl. 



In the upper and lower panels of Figure [U we show the 
confidence contours of the FI and the BI spectra, re- 
spectively. The confidence contours touch at the ^99% 
level of significance for the FI spectra and at the ^^95% 
level of significance for the BI spectra. The probabilities 
that the flux-energy parameters of the first absorption 
line (model 7; Table [3]) are the same between 0BS2 and 
0BS3 (nuU probabilities) are <1 x 10"'' and <2.5 x lO^^ 
in the FI spectra and BI spectra, respectively.^'* To take 
into account possible sampling effects caused by the num- 
ber of trials used in our variability analysis we multiply 
the null probabilities by six. This factor corresponds to 
the number of absorption lines (two) times the number 
of observations (three). We conclude that the variability 
of the first absorption line is significant at the >99.9% 
and >;98% levels in the FI and BI spectra, respectively. 

We note that the slight differences ( < 68% signifi- 
cance) between the FI and BI confidence contours may 
possibly be associated with differences in the responses, 
variations in the signal-to-noise ratios of the two detec- 
tors, and statistical noise. 

3.1.2. XSTAR spectral fits. 

The spectral analysis presented in ^3.1.11 indi- 
cates that the intrinsic X-ray absorbing medium of 
APM 08279+5255 is complex and contains absorbers 
with different properties (see models 5-7 in Table [3|). 
We identified a low-ionization absorber with a col- 
umn density of log A^h'^23 and an ionization parame- 
ter of log ^<0 1^ (see models 2-4 in Table [3]). This 
low-ionization absorber is required to model the absorp- 
tion detected below ~2 keV (observed-frame). An ad- 
ditional complex absorber is required to fit the broad 
absorption features with rest-frame energies between 

epoch OBS2, the first absorption feature (absl) falls near the Si 
K edge where events with energies lying in the range 1.7-1.95 keV 
were ignored in the analysis. Although the loss of these data points 
adds a larger statistical error to the best-fitted parameters of absl, 
it does not significantly affect the analysis. Specifically, this error 
in the best-fitted energies is less than the value of (T^bsi a-nd more 
Ukely close to ~ 0.2 keV (see Table|3J. 

The square of the probabilities of being outside the confidence 
contours that barely touch (see Figure [2]l is an upper limit to the 
null probabilities. 

The spectral fits did not show an improvement using a warm- 
absorber model; however the <4 keV absorption was not well con- 
strained since the Suzaku spectra start at rest-frame energies of 
S ~ 2 keV. 



7-12 keV. Given the range of energies and variability 
of the broad absorption features, we interpret this ab- 
sorption as a blend of highly ionized (2.75<log C<4) iron 
absorption lines blueshifted by an outfiow. This expla- 
nation is consistent with recent model s that attempt to 
simulate X-ray BALs in quasars (e.g.. ISchurch &: Done I 
l2007t l. To test this interpretation we next employ 
more complex, but more realistic, models to fit the 
APM 08279+5255 spectra. 

As a first attempt we fit the low and high-energy 
absorption of APM 08279-^5255 in epoch OBSl with 
a model that includes a power-law (with Galactic ab- 
sorption) and one warm absorber (model XSTARl, Ta- 
ble E]). The warm-absorber model i s calculated using 
the XSTAR code (s ee, e.g.. iKallman k. Bautista 200l|; 
iKallman et"all I1996D . XSTAR calculates the physical 
conditions and absorption-emission spectra of photoion- 
ized gases with variable abundances. In the current anal- 
ysis we use a recent implementation of the XSTAR code 
called WARMABS that can be used as a model within 
XSPEC. For the WARMABS model we assume turbu- 
lent velocities fturb = 1,000 km s""'^ (default velocity of 
the model). For the fits we assumed solar abundances, 
a redshift of 3.91 for the warm absorber, and we left the 
column density and ionization parameter of the warm 
absorber free to vary in the fit. We note that spectral 
fits using model XSTARl attempt to fit both the the low 
and high-energy absorption of APM 08279-1-5255 with a 
single warm absorber. 

To constrain the iron abundance (apc), we allowed this 
parameter to vary in the spectral fits in the model XS- 
TAR2. The only difference between models XSTARl and 
XSTAR2 is that ape is fixed in model XSTARl and free 
to vary in model XSTAR2. We find that allowing apc 
to vary in our spectral fits does not lead to a significant 
improvement in the fits, and the best-fit values of ape 
in fits with model XSTAR2 are consistent with no iron 
over-abundance. Our results do not confi rm an apparent 
iron o ver-abundance (aFc>2) claimed bv iHasinger et ahl 
(j2002j) and Ramirez (2008) based on their analyses of 
previous observations of APM 08279-1-5255. 

We next assumed the spectral model XSTAR3 consist- 

Since a warm absorber wi th log g ^ 1.0 is expected to have 
a temperature of < 10® K (e.g.. IChelouche fc Netzer1[2003) . we do 
not expect a thermal broadening higher than 100 km . 
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TABLE 5 

Results from spectral fits using XSTAR to epochs OBSl, OBS2 and OBS3 of APM 08279+5255. 



XSPEC Model Paramctcr 



OBSl 



OBS2 



OBS3 





Values FI 


Values BI 


Values FI 


Values BI 


Values FI 


Values BI 


r 


1.96±0.04 


1.99±0.05 


2.00±0.05 


1.96±0.07 


1.97±0.05 


1.95±0.06 


log Nil 


23.35±0.08 


23.19±0.10 


23.30±0.09 


23.28±0.12 


23.25±0.09 


23.26±0.08 


log 5 


1.02±0.52 


0.45±0.36 


0.99±0.52 


1.76±0.62 


1.05±0.75 


1.51±0.66 




69.1/71 


70.7/69 


82.0/64 


92.7/64 


73.8/69 


48.2/69 




0.54 


0.42 


0.06 


0.01 


0.32 


0.97 


r 


1.94±0.05 


2.00±0.07 


2.01±0.05 


1.95±0.07 


1.97±0.05 


1.95±0.06 


log A^H 


23.30±0.09 


23.18±0.10 


23.40±0.08 


23.24±0.12 


23.25±0.09 


23.26±0.08 


log 5 


1.09±0.51 


0.54±0.42 


1.28±0.65 


1.70±0.78 


1.05±0.75 


1.51±0.66 


apo 


1.1±0.3 


1.8±0.5 


0.6±0.2 


1.2±0.5 


l.OitO.3 


1.0±0.4 


x^ / y 


68.7/70 


65.2/68 


77.1/63 


92.5/63 


73.8/68 


48.2/68 




0.52 


0.57 


0.11 


0.01 


0.29 


0.97 


r 


1.97±0.05 


1.96±0.08 


2.02±0.05 


1.95±0.07 


1.96±0.05 


1.92±0.06 


log A^Habsl 


23.31±0.09 


23.17±0.18 


23.40±0.08 


23.28±0.10 


23.21±0.09 


23.14±0.09 


log ?absl 


1.17±0.55 


0.50±0.36 


1.26±0.47 


1.28±0.25 


0.78±0.46 


1.31±0.49 




0.8±0.2 


1.5±0.7 


0.6±0.3 


1.2±0.6 


l.litO.3 


1.3±0.5 


Zabsl 


3.91 


3.91 


3.91 


3.91 


3.91 


3.91 


log Arjjabs2 


23.07±0.31 


23.36±0.26 


22.79±0.36 


22.87±0.41 


22.91±0.35 


23.21±0.39 


log ?abs2 


3.7±0.4 


3.6±0.5 


3.5±0.3 


3.8±0.5 


3.5±0.4 


3.7±0.5 


Zabs2 








1 7Q-i_n 1 9 


9 1 Q-i-n 1 ^ 
z. lymu. lo 


9 07-1-0 1 


X / y 


58 0/67 


54 9 /65 


73 6/60 


89 8/60 


70 2/65 


43 6/65 


P(y^ IvY 


0.77 


0.81 


0.11 


8 X 10~^ 


0.31 


0.98 


V 


1.94±0.06 


1.95±0.07 


2.01±0.05 


1.92±0.07 


1.93±0.05 


1.91±0.07 


log A^Habsl 


22.86±0.16 


23.01±0.26 


22.93±0.31 


23.16±0.36 


22.82±0.21 


23.02±0.28 


log ^absl 


3.8±0.3 


3.5±0.2 


3.6±0.2 


3.8±0.5 


3.1±0.3 


3.3±0.3 


^absl 


3.08±0.10 


3.22±0.12 


2.78±0.08 


2.57±0.12 


3.24±0.06 


3.12±0.08 




(0.19±0.02) 


(0.16±0.03) 


(0.27±0.02) 


(0.32±0.03) 


(0.15±0.01) 


(0.18±0.02) 


log AfHabs2 


23.38±0.15 


23.43±0.22 


22.91±0.35 


23.09±0.39 


23.10±0.38 


23.07±0.41 


log Cabs2 


3.6±0.5 


3.4±0.2 


3.4±0.3 


3.7±0.6 


3.5±0.3 


3.6±0.6 


Zabs2 


2.10±0.06 


2.17±0.08 


1.97±0.09 


1.78±0.13 


2.18±0.12 


2.05±0.10 




{0.45±0.02) 


{0.43±0.02) 


(0.48±0.02) 


(0.53±0.03) 


(0.42±0.03) 


(0.46±0.03) 




56.3/66 


55.1/64 


65.1/59 


80.0/59 


67.1/64 


43.4/64 


PixV^Y 


0.80 


0.78 


0.27 


0.04 


0.37 


0.98 



XSTARl. 



XSTAR2. 



XSTAR3.. 



XSTAR4.. 



XSTAR1= XSPEC model warmabs(pow); XSTAR2= XSPEC model warmabs(pow) (FeA variable); XSTAR3= XSPEC model 
warmabs*warmabs(pow); XSTAR4= XSPEC model wabs*zwabs*warmabs*warmabs(pow) ^ ^^^^^ and ^'^^^^ corresponds to the 
estimated outflow velocities. They are calculated from equation [2] based on the redshift of the absorbers of model XSTAR4. 
They are not parameters of the spectral fit.^ Probability that / f is greater than the value obtained. 



ing of a power-law, one stationary ionized absorber with 
a turbulent velocity of wturb = 1, 000 km s~^ and a sec- 
ond outflowing ionized absorber with a turbulent velocity 
of wturb = 10,000 km s^^ (the maximum value allowed 
by the model^^). We allowed the ionization parameters 
of both ionized absorbers and the redshift of the second 
ionized absorber to vary in the fit (model XSTAR3; Ta- 
blets]). We find that the best fitted redshift of the second 
warm absorber is z 2 in both the FI and BI spec- 
tra. The i^-test indicates an improvement in the fits of 
OBSl with model XSTAR3, that assumes two ionized 
absorbers, compared to fits with models XSTARl and 
XSTAR2, that assume a single ionized absorber, at the 
>99.5% of significance level in the FI and BI spectra. 
We conclude that a single warm-absorber model cannot 
accurately fit both the low and high-energy absorption 
in APM 08279+5255. 

^"^ The high-enengy absorption features modeled with Gaussian 
absorption lines in model 7 could be the result of one or more highly 
ionized absorbers. The Doppler broadening velocities of each ab- 
sorption line component of model 7 are ~ g,abBi (where i=l,2 in- 

dicates the component). From Table [3] these Doppler broadening 
velocities are at first order comparable to the assumed values of 

fturb- 



We finally assumed the spectral model XSTAR4 con- 
sisting of an absorbed power-law and two outflowing 
ionized absorbers. We assumed turbulent velocities of 
wturb = 10, 000 km s^"^ for the flrst and second outflow- 
ing ionized absorber (model XSTAR4; Table E]). The 
main difference between models XSTAR3 and XSTAR4 
is that the redshift of the flrst warm-absorber is flxed in 
model XSTAR3 to the systemic redshift of the quasar 
and free to vary in model XSTAR4 and that model XS- 
TAR4 includes a neutral absorber. For fits using model 
XSTAR4, we allow the redshifts, column densities, and 
ionization parameters of the absorbers to vary. The 
values for these fits are similar to those found for model 7 
(see Table [5] and [3]). We find on average best-fit redshifts 

of Zabsl ~ 3 (Uabsl ~ 0.2c) and 2abs2 ^ 2 (Wabsl ~ 0.5c), 

and column densities of log A'^n.absi ^ 23.0 and log 
A'H,abs2 ^ 23.2, where abs\ and ahs2 correspond to the 
two warm absorbers. We confirm the results of model 
7, by finding a significant change in the best-fltted red- 
shift of the flrst warm-absorber component in model XS- 
TAR4 between epochs 0BS2 and 0BS3 (see Table [5]). 
Even though the two warm-absorber model results in ac- 
ceptable flts, the best-fit parameters should only be con- 
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sidered as basic estimates of the wind properties since 
the kinematic and ionization structure of the outflow are 
likely to be more complicated. Our spectral fitting results 
of models that include ionized absorbers (see Table [5]) in- 
dicate that both models XSTAR3 and XSTAR4 provide 
acceptable fits to the spectra of APM 08279-1-5255 for 
epochs OBSl and 0BS3, however, model XSTAR4 pro- 
vides a better fit to the data for epoch 0BS2 than model 
XSTAR3. The i^-test indicates that in epoch 0BS2, 
spectral fits using model XSTAR4 provide an improve- 
ment over fits using model XSTAR3 at the ~ 99% level 
of significance in the Fl and BI spectra. 

3.2. PIN spectral analysis. 

We also examined the spectrum of APM 08279-1-5255 
in the 10-40 keV energy band using the PIN-HXD 
data. Unfortunately, no signal from the source above the 
background in the NXB was found. The background- 
subtracted source spectrum was found to be within 
5% of the non- X-ray backgrou nd (NXB) spectrum pro- 
vided by the Suzaku team (e.g.. lMizuno et aLll2007t l. We 
arrive at a similar conclusion from the count rates pre- 
sented in Table [2l The non-detection with the PIN pro- 
vides an upper limit on the flux density at 20 keV of 
APM 08279-^5255 of - lO-^photons s'^ keV"^ This 
limit is consistent with an extrapolation of the XIS spec- 
trum. 

4. DISCUSSION 

The X-ray spectrum of APM 08279-f 5255 is known to 
contain absorption features at rest-frame energies above 
7 keV (Chartas et al. 2002; Hasingcr et al. 2002). These 
features have been interpreted in the p ast in two dif- 
ferent ways. The first interpretation by IChartas et al.l 
(I2OOI was based primarily on the analysis of the 2002 
Chandra observation of APM 08279-1-5255 and posits 
that the absorption features are due to highly blueshifted 
Fe XXV Ka and/or Fe XXVI Ka absorption fines. An 
XMM-Newton observation of APM 08279-^5255 per- 
formed 1.8 weeks (proper-time) after the Chandra obser- 
vation showed significantly different high-energy absorp- 
tion s tructure which was interpreted by Chartas et a]] 
(|2003t ) to imply variability of the absorption features 
over timesca les of the order of wee ks. The second inter- 
pretation by iHasinger et al.l (|2002f ) is based on the 2002 
XMM-Newton observation of APM 08279-^5255 and pro- 
poses that the high-energy absorption feature arises from 
an iron absorption edge produced by a metal enriched 
(Fe/0 « 2-5 Fe/Oj©) ionized absorber. One impor- 
tant conclusion from section 13.1.21 is that the absorp- 
tion feature found at 7-12 keV rest-frame can be fitted 
with two highly ionized blue-shifted warm absorbers that 
do not require super -solar mctallicit ies. We note that 
IHasinger etall (|2002D and [Ramirez, lj2008f ) had claimed 
iron over-abundance (aFe>2) based on their analyses 
of previous observations of APM 08279-1-5255. In sup- 
port of the two-component "iron-blend-outflow" scenario 
we mention that the two absorption-line model and the 
notch model (models 7 and 5; Table |3|) provide sig- 
nificantly better fits to the absorption feature between 

The background-subtracted source spectrum is obtained by 
subtracting the non-X-ray-background and an estimate of the X- 
ray background from the total detected PIN spectrum 



7-12 keV rest-frame than an absorption-edge model in 
epoch OBSl (see §3). We also note tha t, in the 2002 
~90 ks Chandra observation analyzed in IChartas et al.l 
(2002), the model containing two absorption lines suc- 
cessfully fits the 7-12 keV rest-frame feature whereas an 
absorption-edge model did not provide an acceptable fit. 
The absorption described either by two absorption lines 
or a notch, may crudely represent absorption through 
an outflow with a large velocity gradient along the flow. 
Variability of the kinematic and ionization state of the 
outflow may explain why, depending on the observa- 
tion, these absorption features could be modelled by ei- 
th er a notch, absorp tion lines, or an edge (see figure 4 
of ISchurch fc Done I [20071 ). We conclude that a time- 
variable outflow provides a plausible explanation for all 
the past X-ray obse rvations of the absorption features of 
APM 08279-f 5255 (jChartas et al.1 [20021 : [Hasinger et all 
l2002f ) including those analyzed here. In ^4.1! we provide 
plausible explanations for the observed variability of the 
high-energy absorption. In §4.2! we use the results of 
our spectral analysis to place constraints on the kine- 
matics of the outflow and in ^4.3! we provide estimates 
of the mass-outflow rate and efficiency of the outflow of 
APM 08279-^5255. 

4.1. Origin of the Variability of the High-Energy 
Absorption Feature. 

Assuming the first interpretation of the origin of the 
high-energy absorption features in APM 08279-f 5255 the 
observed shift in the energy of the first absorption line be- 
tween epochs 0BS2 and 0BS3 is likely due to a change 
in the outflow velocity of the absorber. Two alterna- 
tive explanations of the shift are a change in the direc- 
tion of the outflow (with respect to the line of sight) 
and a variation in the ionization parameter of the ab- 
sorber. A change in the direction of the outflow is ex- 
pected to show a shift in energy of both components of 
model 7 (Table [3]). We only find variability in one com- 
ponent (absl, see §3); however this picture could still 
be valid if each outflow component is driven indepen- 
dently. A change in the ionization parameter seems to 
be a less probable scenario. We checked this by fitting 
the spectra of epochs 0BS2 and 0BS3 simultaneously 
with model XSTAR4, keeping the redshift of absl as the 
only common parameter between the fits. In the case 
where a change in the ionization parameter produced the 
detected variability of the first absorption line absl, we 
expect that the fits to the spectra of epochs 0BS2 and 
0BS3 will not be improved by allowing the redshift of 
absl to vary independently in these fits. For the simul- 
taneous fit to the spectra of epochs 0BS2 and 0BS3, 
where the redshift of absl was kept a common parameter, 
we find Zabsi(FI) = 3.11 ± 0.04, Zabsi(BI) = 3.04±0.06, 
xlJiy = 138.1/124 and Xm/^^ = 136.3/124. In the 
case where we fit the spectra of epochs 0BS2 and 0BS3 
independently using model XSTAR4 (see Table [5]) we 
find xli/j^ = 132.2/123 and Xm/^ = 123.4/123.1^ The 
improvement based on the i^-test of fitting the spectra 
of epochs 0BS2 and 0BS3 independently, compared to 
keeping a common redshift of absl, is at the ~99% and 
>99.9% level of significance in the FI and BI spectra. 

These values are obtained by summing the and degrees of 
freedom of epochs OBS2 and OBS3 in Table [5] 
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Fig. 3. — Simulated 6—8 keV absorbed spectra (log(i?FE)) assuming an incident power-law spectrum with a photon index of r=2. 
The ionized absorber is modeled with the XSTAR model WARMABS assuming a column density of log Nn = 23, solar abundances, 
''turb = 1000 km s~^, and the following six values of the ionization parameter log ^=2.75, 3.00, 3.25, 3.50, 3.75, and 4.0, respectively. Some 
of the main absorption lines in this range of ionization states have been marked in the figure. Among them Fe XX Ka (2s^2p^ — ls2s'^2p^; 
6.50 keV), Fe XXV Ka (Is^ - ls2p; 6.70 keV), Fe XXV K/3 (Is^ _ ls3p; 7.89 keV), Fe xxvi Ka (Is - 2p; 6.97 keV), Fe xxvi K/3 (Is - 3p; 
8.27 keV). 

and i?max — £'abs2 + 2(Tabs2- From the best-fit val- 
ues of the APL-fNo model (model 5; Table [3]) we have 

£^inin = £^notch — W^notch/2 and i?max = ^^notch + Wnotch/2. 

The values of Emin and i^max are presented in Tabled In 
this table the values of i?min and i?max are shown sepa- 
rately for the FI and BI spectra and for the two different 
models used to obtain them (model 7=APL-|-2AL; model 
5 =APL+No). 

The relative strengths of the iron resonance absorption 
lines will depend on the ionization state of the outflow- 
ing medium. To demonstrate this effect in a basic way 
we have performed several simulations using the warm- 
absorber model (WARMABS) of XSTAR. In Figure [3] 
we show a simulated absorbed spectrum in the 6-8 keV 
rest-frame energy range for an absorbing medium with 
solar composition, log iVH=23, and having six different 
values of the ionization parameter (2.75 < log ^ < 4.0). 
The two strongest iron lines for this highly ionized ab- 
sorbing medium have rest (or laboratory) energies of 6.70 
keV (Fe xxv Ka; Is^ - ls2p) and 6.97 keV (Fe xxvi 
Ka Is — 2p). In general the Fe xxv Ka line will be 
stronger than the Fe xxvi Ka line for a medium with 
2.75 < log ^ < 4.0. Therefore the absorption at the lower 
end of the 7-12 keV rest-frame range is most likely as- 
sociated with the Fe xxv Ka line. Based on our inter- 
pretation of the high-energy absorption features we esti- 
mate Vjnin assuming that the absorption at the low end 
of the X-ray BAL is due to a line arising from highly 
blueshifted Fe xxv Ka (Ei^^, = 6.7 keV). On the other 
hand to estimate Umax we make the conservative assump- 
tion that the absorption at the high end of the X-ray BAL 
is due to a line arising from highly blueshifted Fe xxvi 
Ka {Eiah — 6.97 keV). The estimated values of Vmin and 
t^max obtained through the procedure outlined above are 



We conclude that the variability of the energy of the 
first absorption line absl is likely not driven by changes 
in the ionization parameter of absl. Suzaku cannot re- 
solve the images of APM 08279-1-5255, however, the time- 
delays between the two brightest images A and B of 
APM 08279-F52 55 is estimated to be of the order of a 
few hours (e.g.. iMun oz et aT1l200 1f). much shorter than 
the observed variability of the high-energy absorption 
feature. We therefore do not expect the combined X- 
ray spectrum of all images of APM 08279-1-5255 to differ 
from that of the individual images within the time-delay. 

4.2. Constraints on the Kinematics of the Outflow. 

Under the premise of the outfiow interpretation to ex- 
plain the absorption at rest-frame 7-12 keV we expect a 
continuous distribution of outflow velocities. This range 
of velocities leads to the the Doppler shift of the energies 
of the resonance absorption lines. The absorption-line 
rest-frame energies Ei^h will thus be shifted to the ob- 
served energies Sobs according to 



Eiab/Eobs = 7(1 - (3cos6), 



(2) 



where 7 is the Lorentz factor, 9 is the angle between the 
wind and our line of sight (l.o.s), and (3 — v/c. 

The minimum and maximum projected velocities 
{v-aam 1'max) of thc outfiow are estimated from the mini- 
mum and maximum energy ranges (i^min, £^max) of the 
high-energy absorption features in APM 08279-1-5255. 
We obtained E-axm and i^max from our spectral fits assum- 
ing first the two absorption-line (APL+2AL) model and 
second assuming the notch (APL+No) model. Specifi- 
cally, based on the best-fit values of thc APL+2AL model 
(model 7; Table [S]), we obtain i?min = £^absi — 2f7absi 
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TABLE 6 

The minimum and maximum energies and velocities of the high-energy 
absorption features in apm 08279+5255. 



Model" 


OBS 


Instrument 


-^min 


-E'max 


^min 










[keV] 


[keV] 


[c] 


[c] 



5 


1 


XIS FI 


7.60±0.24 


11.58±0.24 


0.13±0.03 


0.52±0.02 


7 


1 


XIS FI 


7.52±0.36 


12.16±0.52 


0.12±0.05 


0.56±0.04 


5 


2 


XIS FI 


7.82±0.37 


11.68±0.37 


0.16±0.05 


0.52±0.03 


7 


2 


XIS FI 


7.88±0.48 


11.57±0.51 


0.17±0.06 


0.51±0.04 


5 


3 


XIS FI 


7.38±0.34 


11.42±0.34 


0.10±0.05 


0.50±0.03 


7 


3 


XIS FI 


7.51±0.32 


11.70±1.06 


0.12±0.05 


0.52±0.08 


5 


1 


XIS BI 


7.46±0.53 


12.17±0.53 


0.11±0.08 


0.56±0.04 


7 


1 


XIS BI 


7.24±0.51 


12.47±1.12 


0.08±0.07 


0.58±0.08 


5 


2 


XIS BI 


7.91±0.78 


12.92±0.78 


0.18±0.10 


0.61±0.05 


7 


2 


XIS BI 


8.74±0.58 


12.24±1.16 


0.28±0.07 


0.56±0.08 


5 


3 


XIS BI 


7.70±0.63 


12.22±0.63 


0.15±0.09 


0.56±0.04 


7 


3 


XIS BI 


7.53±0.82 


12.40±0.92 


0.12±0.11 


0.57±0.06 


7 


Cha02'' 


ACIS BI 


7.95±0.11 


10.28±0.22 


0.18±0.01 


0.40±0.02 


7 


Has02'' 


EPIC pn 


6.95±0.44 


15.00±1.06 


<0.10 


0.72±0.05 



^ Model used to estimate Emin and Emax- Model 5 (see Table [sj is a power-law with Galactic 
absorption, intrinsic absorption, and a notch absorber; Model 7 (see Table [3} is a power-law 
with Galactic absorption, intrinsic absorption, and two absorption lines. ^ In this Table we 
identify as Cha02 the 88.8 ks observation of APM 08279-1-5255 performed with Chandra in 
2002 and analyzed in detail in lChartas et al.| (120021) . We also identify as Has02 the 100.2 ks 
observation of APM 08279+5255 performed with XMM-Newton in 2002 and analyzed in detail 
in IHasinger et al.l II2002D . 
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Fig. 4. — Wind velocity (in units of Qjc) plotted as a function of radius from the central source (left panel) and as a function of time 



(right panel) for a radiation-pressure driven wind. Notice that the radius is in units of the Schwarzschild radius [Rs 



-) and the time 



is in units of - 



The wind velocities are calculated for _Ri, 



: 3-Rs and Ri- 



IbRs, respectively. 



presented in Table [6l These velocities are obtained us- 
ing equation ([2|) assuming that our line-of-sight makes 
an angle of = 20° with the velocity of the outflow 
(jChartas et al.ll2007af) . A change of 10° in 9 will intro- 
duce a variation of <10% in our estimates of the veloc- 
ities. The velocities of the outflow range between G.lc 
and 0.6c. The mean value of i?max from all CCDs and 
observations is 12.04±0.22 keV corresponding to a mean 
value of Wmax = 0.55±0.02c. This maximum value of the 
outflow velocity constrains the angle between the outflow 



direction and our hue of sight to be < 36°. This rela- 
tively small angle is consistent with outflow models that 
posit that BAL quasars are viewed through coUimatcd 
outflows. 

As argued in lChartas et al.l (|2002ll2003l . l2007aD we are 

likely observing the X-ray absorbers as they are accel- 

The Doppler-shift formula (equation [2]l predicts that given 
a fixed ratio of i?iab/£^obs=^lo the maximum angle between 
our line of sight and the wind direction is given by Smax = 
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erated near their launching radii. We use the foUowing 
equation to describe the basic dy namics of a radiatiq n- 
driven outflow (e.g., equation 1 of lChartas et al.ll2002f ): 



.d = [(Tf 



Luv 



-!)(■ 



1 



1/2 



(3) 



^Edd ^launch 

where Wwind is the outflow velocity in units of c, Tf is 
the force multiplier, LEdd is the Eddington luminosity, 
^launch is the ladius (units of Rs) at which the wind 
is launched from the disk, and R is the distance (units 
of Rs) from the central source. The expression for the 
dynamics of the outflow can be simplified by defining 

^/ = - 1. 

In Figure 2] we plot wind velocity versus radius (left 
panel) and wind velocity versus time (right panel) for an 
outflow launched at radii of 3Rs (solid lines) and 15i?s 
(dashed lines). In Figure|3]the units of velocity, distance, 
and time are r2/c, Rs, and Rs /(^f c), respectively. As- 
suming Luv 0.2Lboi (Irwin et al. 1998' ), Lboi/LEdd 
0.2 a nd T/ - 100 (e.g.. Arav et al. 1994 : Laor fc Brandt I 
|2002[) . we have Vlf ~ 1.7. Lboi/iEdd ~ 0.2 is ob- 
tained using the Lb oi/I/Edd vs. F c o rrelation found 
in RQ quas ars (e.g., IWang et al.1 120041 : IShemmer et al.1 

[200l [200l fo r F ~ 1.9 T herefore with Lboi = 

TxlQiVZ^^O Ulrwin et al.lll998l : lReichers et al S^MS)^^ 
Mbh - IQI^Mq^^^ where hl ~ 100 (jEgami a l. 2000)^2 
is the lens magnification factor. Equation [3] could be 
modified using a reliable SED describing the central 
source, adding relativistic corrections, and calculating 
the force multiplier at every point of the trajectory of 
the outflow. We stress, however, that our simplified ap- 
proach is sufficient to provide first order approximations 
to the launching radius and the time scales involved in 
the dynamics of the outflow. Equation [3] can be writ- 
ten as i?iaunch/-Rs = f^/(c/woo)^; therefore for Qf ^ 1.7 
-Riaunch ~ 3 X {c/vooYRs- The latter expression can be 
used to obtain first order approximations of the launch- 
ing radius given the velocity of the outflow. 

Assuming a radiation-driven wind, it is expected that 
the time required to accelerate the outflow to fractions 
of c is of the order of 10-^ (see Figure |4]). For the 
black-hole mass of APM 08279-f 5255 of A/bh ~ W^^Mq 
we estimate that the time to accelerate an absorber to 
near-relativistic velocities is ~weeks (rest-frame). We 
have reported in this work probable variability of the 
high-energy absorption features over a time-scale of ~1 
month (rest-frame). This short time-scale variability is 
consistent with the expected variability timescale of a 
radiation-driven wind. 

4.3. Constraints on Mass- Outflow Rate and Efficiency 
of the Outflow. 

Based on our estimated values of the outflow veloc- 
ities, column densities, and launching radii we present 
constraints on the mass-outflow rates and outflow effi- 
ciency associated with the outflowing X-ray absorbers of 

We note that an estimation of Lbol based on the optical and 
UV spectra should be more precise than a n extrapolation of Lbol 
based on X-ray luminosities as it is done in lRami'red ()2008l l . 

See, however, Reichors ct al. ( 2003) that find a magnification 
of fj,]^ 4. iReichers et al., (,2008 ) also use the observed width of 
the CIV line to obtain a black-hole mass of Mbh ~ ^^^^ l^i, 



APM 08279-f 5255. The efliciency is deflned as the ratio 
of the rate of kinetic energy injected into the ISM and 
IGM by the outflow to the quasar's bolometric luminos- 
ity, i.e., 

ex = i-T— where M = inR^pyfc = i-Kfc^-^NumpV, 

(4) 

where fc is the covering fraction, A'h is the column den- 
sity, R is the radius, and Ai? is the thickness of the 
absorber. To estimate the efficiency we use the two 
absorption-line model (APL-I-2AL; model 7 of Tabled. 
We calculate the bulk velocities of each outflow com- 
ponent based on the energies of the absorption lines 
in model 7 and through the use of equation (|2 l) with 
_^ab=6.7 keV and 6 = 20°. As in IChartas erair(|2002l 
l2003f ) we interpret the high-energy absorption features 
as being due to highly ionized Fe (Fe xxv Ka) in a gas 
with solar abundances, and we estimate log A'h using a 
curve-of-growth analysis. In Table [7] we present the out- 
flow velocities, Wabs, the column densities, log N^i, the 
mass-outflow rates, M, and the outflow efficiencies, ck, 
of the two modeled absorbers of the outflow. We note 
that the values of the column densities and velocities in 
Table [7] are consistent with those found using the pho- 
toionization code XSTAR (see Table [5]). 

To obtain error bars for tK and M we performed a 
Monte Carlo simulation, assuming a uniform distribu- 
tion of the parameters fc, R and R/AR around the ex- 
pected values of these parameters, and a normal distri- 
bution for log iVn (described by the parameters in Ta- 
ble [7]) . Specifically, we assume a covering factor lying in 
the range fc — 0.1 — 0.3, based on the obse rved fraction of 
BAL quasars (e.g.. lHewett fc Foltz |[2003l ) and a fraction 
R/AR ranging from 1 to 10 bas ed on current theo retical 
models of quasar outflows (e.g. JProga et al]|2000r ). Note 
that in Table [7] we also include the outflow parameters 
of the -90 ks Chandra observation of APM 08279+5255 
performed in 2002 (|Chartas et al.ll2002[ ). Based on our 
estimated maximum outflow velocities (wmax ~ 0.6c) we 
expect that R will be similar to -Riaunch and range be- 
tween 3Rs and 15i?5 (see Figure H]). We note that this 
is a conservative assumption since larger values of R 
will result in larger mass-outflow rates and larger effi- 
ciencies. Additionally, the short variability time-scales 
(—weeks) are also consistent with a launching radius of 
a few times Rs- Variabihty in APM 08279+5255 over 
time-scales —weeks has been previously reported (Char- 
tas et al. 2003) based on the differences of the high- 
energy absorption features detected in the Chandra and 
XMM-Newton observations. Our current analysis of the 
Suzaku observations of APM 08279+5255 allows us to 
compare absorption features observed with the same in- 
struments; therefore, it avoids any possible systematic 
uncertainties due to differences in the responses of the 
instruments. 

Our results indicate that the average fraction of the to- 
tal bolometric luminosity of APM 08279+5255 injected 
into the IGM in the form of kinetic energy is ex = 0.7 ± 
0.3. From the results of Table[7]we also obtain the aver- 

The velocities obtained in Table [5l assume the redshifts of the 
absorbers in model XSTAR4 are due to the relativistic Doppler 
effect (see equation [2]l . 
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TABLE 7 

Projected maximum outflow velocities, mass-outflow rates and efficiencies of outflows in APM 08279+5255 



OBS 


Instr. 




log AfH(absl) M (absl) (absl) 




log 7VH(abs2) M (abs2) (abs2) 












[Mef^l^yi-^] 



1 


XIS FI 




22.84±0.19 


470+478 
4'U_3ig 


'-'■'J^-o.oi 


n 47+0-02 
'^■^'-0.03 


23.07±0.27 


9079 "I" ^"^^^ 
^■^'^-1496 


5+°-^ 
U.0_o.4 


2 


XIS FI 


97+0.03 
-0.02 


22.86±0.25 


71 q+794 
'^■^-507 


n 05+° "'* 


n 40+0.02 


22.81±0.39 


14681?;^! 




3 


XIS FI 


-0.01 


22.87±0.28 




'-'■'Ji-0.01 


n 44+0-04 


23.05±0.19 


1678llI°I 


U.0_o 2 


1 


XIS BI 


Q ,7 + 0.03 


22.89±0.24 


475+^3° 


'J-'Ji-o.oi 


n 46+°°3 
'^•^O-o.oa 


23.48±0.20 


4772lf,i^ 




2 


XIS BI 


34+0.02 


23.01±0.36 




iq+0-25 
'-'■^^-0.15 


n ,ro + 0.05 


23.09±0.36 


2849l382g 


o.8iJ:^ 


3 


XIS BI 




23.10±0.21 


gg2 + 1006 

°" -663 


U-UO_o 03 


n "sn+''-°'' 

'^•"'^-0.04 


23.14±0.24 


24861^722 


7+°'^ 

"■'-O.B 


Ch02 


ACIS BI 




22.99±0.11 


g22+570 


'J-lJ'3_o.02 




23.07±0.12 


1507t«f 




Has02 


EPIC pn 


21+0 03 


23.10±0.36 


11^:^+1550 


10+°" 


n 50+0-06 


23.36±0.30 


450015443 


1-^-n.c) 



^ The estimated values of the outflo-w properties were based on fits 
The values of M and ek are obtained by equation Q assuming Mbh 
IReichers et al.|[200l) . 

age mass-outflow rate of M = 3324 ±915 /x^^Moyr"^. 
On the other hand the mass-accretion rate is Mace = 
^ 484O/Lt^"^M0yr~^ (assuming 77 w 0.1 and Lboi = 

7 X IO^^/^^^Lq). Therefore the mass-outflow rate is 
comparable to the accretion rate. In the context of re- 
cent i no dels of structure f ormation fe.g.. iGranato et all 
[2OOI ISpringel et al.l[200l . our estimated values of e/^ 
in APM 08279-1-5255 suggest that these outflows should 
be an important source of feedback in their host galaxies 
and also play an important role in regulating the growth 
of the central black hole. 

To obtain an independent estimate of the mass-outflow 
rate M, we derive nR^ based on the deflnition of ^ = 
L/{nR^). Assuming log ^ - 3.5 (e.g., model XSTAR4, 
Table [5]) and an ionizing luminosity similar to the X-ray 
luminosity Lx 4 x lO^'^/i^^erg s~^, we obtain nR^ 
1 X 10"'^/i^^ cm^^. Therefore assuming an overall velocity 
of the outflow ~ 0.45c we find M ~ 9 X lOVZ^^oyr"^- 
The value of M derived from the best-fit ionization pa- 
rameter is comparable (within a factor of three) to the 
value found for the second absorber abs2 in Table [71 
where we have estimated the location of the absorber 
from variability arguments. We caution, however, that 
the estimation of the mass-outfiow rate from the ioniza- 
tion parameter assumes a spherical outfiow illuminated 
by a point source. However, since we expect that the 
X-ray absorber is located a few Tg from the X-ray source 
the point source approximation may not be accurate in 
this case. 

5. CONCLUSIONS 

Our analysis of three long Suzaku observations of 
the BAL quasar APM 08279-1-5255 indicates strong and 
broad absorption at rest-frame energies of <2 keV (low- 
energy) and 7-12 keV (high-energy). Based on the 
J^-test the low-energy absorption is significant at the 
>99% and >99.9% levels in the front-illuminated (FI) 
and back-illuminated (BI) Suzaku XIS spectra, respec- 
tively. The high-energy absorption is significant at 
>99.9% (FI spectrum) and at >99% (BI spectrum) con- 
fidence, respectively. The medium producing the low- 
energy absorption is a nearly neutral absorber with a 
column density log A'h'^23. The medium producing the 
high-energy absorption appears to be outfiowing from 
the central source at near-relativistic velocities and with 



that assumed an absorbed power-law model with two absorption lines. 
~ lO^^/x^^M0 (see 31:21 and Ltol = 7 X W^'^ Lq l llrwin et aT]|1998l : 



large ionization parameters (2.75 < log ^ < 4.0), consis- 
tent with results obtained from a previous Chandra ob- 
servation of this object (Chartas et al. 2002). Simula- 
tio ns of highly ionized near -relativistic winds performed 
by ISchurch fc Done I ()2007l ) indicate that the resulting 
X-ray broad absorption profile may have the apparent 
shape of an absorption edge, a notch, or a combination 
of absorption lines depending on the assumed dynamics 
and degree of ionization of the outflowing absorbers. 

Our observations of the 7-12 keV rest-frame features 
are well described by a two component absorber model. 
We find that in epoch OBSl, spectral fits with the two 
component absorption-line model (APL -I- 2AL) are sig- 
nificantly better (>95% confidence level) than fits with 
absorption-edge models (APL + Ed) . We note that spec- 
tral fits with models that included ionized absorbers with 
free iron abundances (models XSTAR2 and XSTAR3 in 
Table [S|) are consistent with no iron over-abundance in 
aU SuzaJfu observations of APM 08279-^5255. 

Our interpretation, of a near-relativistic outflowing ab- 
sorbing medium in a high ionization state (2.75 < log ^ < 
4), is consistent with our analysis of all the past X-ray 
observations of APM 08279-1-5255. Our spectral anal- 
ysis indicates that the outflow velocities of the highly 
ionized absorbers detected in the three Suzaku observa- 
tions range between 0.1c and 0.6c. The maximum de- 
tected projected outflow velocity of ~0.6c constrains the 
angle between our line of sight and the wind direction 
to be ^36°. We find possible variability of the high- 
energy absorption lines between epochs 0BS2 and 0BS3 
at the ^99.9% and ^^98% significance levels in the FI and 
BI spectra, respectively. Our spectral analysis indicates 
that the variability is likely due to a change in the outflow 
velocity of the absorber. The short time-scale (~month 
in the rest-frame) of this variability is probably indicat- 
ing that this absorber is strongly accelerated. This short 
time-scale variability combined with the high ionization 
of the absorbing material imply that the absorbers are 
launched from distances "^IQRs from the central source. 

Assuming our interpretation that the absorption lines 
detected at rest-frame energies of 7-12 keV are due 
to Fe XXV, we estimate that a signiflcant fraction 
(0.7±0.3) of the total bolometric energy over the quasar's 
lifetime is injected into the intergalactic medium of 
APM 08279-1-5255 in the form of kinetic energy with a 
mass-outflow rate of M = 3324 ± 915 ^^^Moyr~^ 
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